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ABSTRACT 

Aims. The aim is to study the stellar content of the bar region to constrain its formation and evolution. 

Methods. Line-strength indices in the bar region of a sample of 6 barred galaxies were employed to derive age and metallicity 
gradients along the bars using stellar population models. 

Results. We find clear radial gradients in the line-strength indices for all the galaxies. We find positive gradients within the bar region 
in the metal indices in four of the six galaxies, and opposite trends in the other two. These latter two galaxies are classified as SAB 
and they present exponential bar light profiles. For all the galaxies we find a positive gradient in the Balmer indices. There is a clear 
correlation between the position of morphological features inside the bar region with changes in the slope and value of the indices, 
which indicate, using stellar population analysis, changes in the stellar populations. Therefore, it seems that the bar regions show a 
gradient in both age and metallicity, changing radially to younger and more metal rich populations for all the galaxies except for those 
two with exponential profiles. 
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,H , 1. Introduction ing dBinnev & Tremainel 1 19871) . so one w ould not expe ct the 

cZ) same abundance trends in both components. Friedli ( 1998|) pre- 

■ The importance of bars in galaxies as mechanisms to transfer dictedi ufjing N . body simulations of bars with pre -existing ex- 

> ; angular momentum and matter on large scaleswithin the galaxy ponential abundances, a null evolution of the stellar abundance 

. r-j is indisputable dKormendy & Kennicuttt |2004l and references profile (although with a deC rease in the mean metallicity) while 

^ therein). Bars are thought to be formed through spontaneous disk the gas abundance profile flatte ned rapidly. As it has been men- 

H ; instabilities or by instabilities produced during galaxy encoun- tioned beforCi the gas . phase abu ndance distribution predicted by 

. , ters; the details of bar formation and their subsequent evolution FriedH , 1998) has been shown to agree with the obsel - V ational 

remains a matter of debate. Until now, only detailed analyses data of abundances in HII reg ions. As for the stellar compo- 

of the gas-phase abundances along bars have been undertaken, nent> there has not been any smdy addressing the stellar pop . 

from observations of emission lines of HII regions. Studies of u i a ti on along bars u sing spectroscopic information. Recently, 

the gas-phase abundances provide us with a present-day snap- Moorthy & Ho ltzman (2006) have presented results for the stel- 

shot of the interstellar medium abundance. These studies have lar popu i ations of bu i ges of spira i ga i axie s, some of which are 

shown that there is little variation in the chemicaljtondances balTed galaxies , finding a difference between bulges of barred 

in the gas-phase along the bar dMartin & FriedhL 119971 119991) . and un . barred ga laxies 
The study of the stellar ages and metallicities, on the other hand, 

gives us 'archaeological' clues as to the formation and evolution To date, it has been difficult to obtain radial abundance dis- 
of the bar. Furthermore, the evolution of the two components, tributions along bars because, in order to obtain reliable stellar 
gas and stars, suffers from very different evolutionary processes; abundances and ages, very high signal-to-noise (S/N) spectra are 
the gas is mainly dominated by the gravitational torque of the necessary, and although bars are high surface brightness struc- 
non-axisymmetric mass component, while the evolution of the tures, this still implies long integration times on medium-size 
stellar component is mainly affected by different orbital mix- telescopes. We have carried out a project to obtain high S/N data 
along the bar major axis of a sample of barred galaxies. This is 

* Based on observations obtained at Siding Spring Observatory the first detailed study of stellar line indices along the bar re- 
(RSAA, ANU, Australia) gion. In this letter, we show clear evidence for line-strength in- 
** Veni Fellow dex gradients along the bars and a correlation between the posi- 

Associate Researcher tion of morphological features inside the bar region (bulge/inner- 

+ Marie Curie disk/inner-bar/primary-bar) and breaks in their respective stellar 

* Retorno JA Fellow populations. 
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2. Observations and data reduction 

2.1. Sample selection 

We have selected barred galaxies from the Third Refe rence 
Catalogue of bright galaxies (RC3) (Ide Vaucouleursl [T948). with 
the following criteria; to be southern (<5 < 0°), classified as 
barred and with inclinations between 10° and 50°, and nearby 
(cz < 2500 km s _1 ) to be able to properly resolve the bar, with 
morphological types earlier than SBb to avoid morphological 
dependency of the results (in this first approach) and to ensure 
that the requested S/N is reached, since early-type galaxies have 
higher bar surface brightness, crucial for the line-strength deter- 
mination. We finally observed six galaxies (names and morpho- 
logical types are indicated in Fig. [TJ. It is by no means a statis- 
tically complete sample, but it is a starting point in the study of 
line-strength indices in bars. 

2.2. Observations and data reduction 

We obtained long-slit spectra of six barred galaxies with the 
Double Beam Spectrograph (DBS) on the 2.3m telescope at 
Siding Spring Observatory (SSO) during February 2006. SITe 
1752 x 532 CCD detectors were used. The gratings employed 
were the 600B and the 600R for the blue and red arms, respec- 
tively, with a slit-width of 2 arcsec. This set-up gives a dispersion 
of 1 . 1 A/pixel for the blue arm and 1 .09 A/pixel for the red arm in 
the wavelength intervals 3892-58 14 A and 5390-73 14 A respec- 
tively, giving a velocity dispersion resolution Act ss 130km s _1 . 
The slit was placed along the bar major axis. The slit length is 
7 arcmin, which allows one to study nearby galaxies in detail, 
while still providing sufficient sky coverage. The position an- 
gle was derived using the DSS images of the sample galaxies. 
Comparison arc lamp exposures were obtained for wavelength 
calibration. Spectrophotometric standards were observed with a 
slit width of 6 arcsec. Additionally, we observed 1 1 G-K stars 
from the Lick/IDS library to be used as templates for velocity 
dispersion (<x) measurements, as well as to transform our line- 
strength indices to the Lick system. The total integration time 
for each of the galaxies was typically 3 hours. 

The data presented in this letter refers to the blue range alone. 
All the spectra were reduced using standard IRAB3 routines. 
Overscan and bias were subtracted from all the frames. No dark 
subtraction was done due to the low dark current of the chips 
used. Flatfielding correction was achieved to the 2% level. The 
relatively small, compared to the slit length, size of the galax- 
ies allowed for accurate sky subtraction for each frame prior to 
combining the reduced spectra. A first order polynomial was fit 
along the spatial direction and then subtracted from the frames. 
Spectrophotometric standards were used to relatively flux cali- 
brate the spectra. For each fully reduced galaxy frame, a final 
frame was created by extracting spectra along the slit, binning 
in the spatial direction to guarantee a minimum signal-to-noise 
ratio of 20 per A in the spectral region of Mgb. 

3. Line-strength indices 

In order to analyse the stellar popu lation of the bars, spec- 
tral indices of the Lick/IDS system (Fab er et ail [l985) were 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agreement 
with the National Science Foundation. 



measured and are presented in F ig. [JJ We adopted the defi- 
nition of iTrager etail d!998l) and IWorthev & Ottavianil d 19971) 
for the high-order Balmer lines H6 and FTy. The Fe3n index 
derived here is a combination of three prominent Fe lines ( 
Fe 4383+Fe5270+Fe5 335)/3 ), similar to the Fe3 index defined 
bv lKuntschnen(l2000h . with Fe4383 instead of Fe5015. 

To compare between different regions of the galaxy and to 
compare with other galaxies, indices have to be measured at 
the same resolution. As to allow flexibility in the model selec- 
tion and to potentially compare with the measurements of other 
galaxy components from other studies, we decided to follow this 
approach. After broadening to the Lick resolutio n, line-strengt h 
indices were measured using the INDEX routine (Cardiel, 1999). 

To derive errors in the line-strength ind ices, we followed a 
similar approach as in Kels on et alj j2006). This method con- 
sists in calculating the errors in the mean flux within a passband 
using noise spectra derived by subtracting a best-fit model to 
the galaxy spectra. To d erive the best-fit we used the OPTEMA 
algorithm described bv lGonzalezI d!993l) .We used the synthetic 
spectral energy distributions provided by Vazdekis et al. (2007). 
Once the error in the m ean flux of each band pass is calculated, 
the formulae derived in ICardiel et alj (1998) were used to cal- 
culate the errors in the indices. As the values of the indices de- 
pends on the Doppler broadening, the velocity dispersion along 
the radii was also calculated. To derive cr and radial veloci- 
tie s, we use t he MO VEL and OPTEMA algorithms described 
by iGonzalezi ([1993). To build an optimal template we use stars 
observed with exactly the same instrumental configuration as the 
galaxies. To correct the indices for the velocity dispersion broad- 
ening we followed the standard procedure of deriving correction 
curves from artificially broadened spectra. We derive, for each 
galaxy, a different curve using the best-fit spectra obtained in 
the derivation of the errors. The reason to derive a polynomial 
for each galaxy is that th e broadening corre ction depends on the 
strength of the indices (Kun tschnen, Il998). and by adopting a 
single polynomial for all the gala xies, artificial tren ds between 
indices and cr can be introduced dKelson et al.Ll2006h . 

To correct the ind ices for any nebu lar emission contribution 
we used GANDALF (ISarzi et al.U 2006) which is a simultaneous 
emission and absorption lines fitting algorithm. We have marked 
all the points with detected emission (see Fig. [TJ, where the 
detection threshold has been calculated with the prescriptions 
given in ISarzi et al.l (12006b . The sensitivity vary from 0.3 A in 
the internal parts to 0.6 A in the external parts of the galaxy. 
When the emission completely fills all the absorption Balmer 
lines, the correction is more uncertain. For those cases, we have 
used different symbols in Fig rjj 

We used TMB03. lBruzual & Charlol d2003l) and Vazdekis et 
al. (2007) models to derive simple stellar population ages and 
metallicities using the indices H/3 and [MgbFe]' (TMB03). We 
also derived [a/Fe] using Fe4383 and Mgb by fixing the ages 
obtained in the [MgbF]' vs. H/3 diagram and the TMB03 models. 



4. Location of the morphological features 

In order to locate the radius of changes in morphological struc- 
ture, ellipse fitting to the light distribution was performed using 
the IRAF task ELLIPSE. All the galaxies showed point-like nu- 
clear regions, so the centre was fixed using the coordinates ob- 
tained by fitting a Gaussian to the nucleus. The position angle 
(PA) and the ellipticity (e) were left as free parameters in the 
fitting. Changes in both ellipticity and PA were considered as 
candidate structures, then visual inspection of the images was 
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performed to see if ther e was any ano maly. The analysis is sim- 
ilar to that described in Erwiji| d2004l) . The images used for the 
analysis come from different sources, all the analysis was done 
with 2MASS data and also a different band in the cases were 
there were published images. The following sources were used 
for each of the galaxies: for NGC 30 81 the HST image o f the 
nuclear region in the F606W band (Malkan et al., 1998), for 
NGC 1433 the fi-band bv lHameed & Devereuxl (1 19991) . and for 
NGC 5101 fl-band image by lEskridge et all dl999l) . We com- 
pared the results in the different bands finding good agreement 
in the position of the morphological features. 

For three out of the six galaxies (NGC 1433, NGC 3081, and 
NGC 2217) there are published data on the radii of m orpholog- 
ical structures within the inner region dErwinl |2004). There is 
good agreement (di fferences of ^ 5%) between our analysis and 
the values found by lErwinl d2004h . 

5. Results and summary 

Figure Q] shows the radial trends of two indices for each of the 
galaxies, a metal and a Balmer index. The other Balmer in- 
dices follow the same trends shown in Fig.Q] The Mg indices in 
some cases differ from the Fe indices, following opposite trends. 
However, one has to bear in mind that raw line indices are af- 
fected by the age-metallicity degeneracy. Therefore, to disen- 
tangle age and metallicity effects from the indices, one requires 
the use of stellar population models. In the cases where H/J is 
more affected by emission H6 was chosen, since it is less affected 
by the emission lines (see Section 3 about emission correction). 
One should be careful deriving the Balmer indices affected with 
emission since the derived values are sensitive to the emission 
correction (increasing the value after the correction). We have 
tried several methods (even no correction) and the trends do not 
change. In Fig. Q]one can clearly see a different trend of the ra- 
dial distribution of the indices with a change at the radius of the 
inner structure (in Fig.Q] the most inner vertical line). 

For clarity, we divide the results of the radial distribution in: 
a) the inner structure, up to the most inner vertical line in Fig.Q] 
and b) the bar region, between the two vertical lines in Fig.Q] 

a) Inner region. There seems to be a gradient of the metal 
and Balmer lines within the inner structure, without taking into 
account the central points. These very central points are possibly 
dominated by the bulge, and show a different behaviour from the 
bar and the inner structure. Both sides of the bar show similar 
trends, reinforcing the results. 

b) Bar region. Within the bar region the slope changes for 
both the metal and the Balmer indices. For NGC 5101 and 
NGC 4643, there seems to be a point where the slope changes 
again corresponding to the radius around the point of maximum 
ellipticity, which is not considered as the end of the bar in the 
ellipticity-PA analysis. NGC 1433 shows an Fe3n distribution 
also compatible with no radial change in the Fe3n, but there is 
a clear difference between the inner part and the bar region. The 
Balmer indices show a positive gradient and a change compared 
to the inner region. For two cases (NGC 2665 and NCG 3081) 
there seems to be a negative trend of the metal indices and a 
slightly negative trend in the Balmer indices; they do not show a 
change in the slope in the bar region; interestingly enough these 
are the only two galaxies classified as SAB in the sample, and 
they present an exponential profile in the bar, which might indi- 
cate a different for mation process and history compared to thos e 
with a flat profile dCombes & Elmegreenill993tlNoguchii ri996). 
This point deserves further investigation. For all galaxies, there 
is also good agreement between both sides of the bar. 



There is an obvious influence of the different components 
to the shape of the radial distribution. We have made tests sub- 
tracting a bulge component (from the decomposition of the light 
profile), this only helps to obtain an even clearer result for the bar 
and inner structure, without changing the basic trends. We prefer 
to show the complete radial profile in order to avoid affecting the 
the observed profiles with any model-dependent analysis. 

Independently of the selected models (see Section |3), we 
derive positive age and metallicity gradients in the bars (outer 
parts of the bars being younger and more metal rich than the in- 
ternal parts), except for the two galaxies with exponential bars 
(NGC 2665 and NGC 3081). In a future paper we plan to do a 
more detailed analysis of the star formation history and chemical 
evolution of all the objects. There, we will present the compari- 
son of these results with dynamical and stellar population mod- 
els using the entire spectra to try to interpret the results in the 
context of bar formation and evolution. 

Summarising, we have obtained high S/N long-slit spectra 
along the bars of six early-type barred galaxies and we have de- 
rived radial profiles of the line-strength indices for the sample 
galaxies. There is a clear difference between the indices in the 
bar region (in Fig. Q] the region between the two vertical lines) 
and the central component. Within the bar region there is a clear 
gradient in both metal and Balmer lines that when combined 
with stellar population models for all galaxies gives a gradient 
in both, metallicity and age, indicating that outer parts of the bar 
are younger and more metal rich. These results might indicate a 
secular inside-out growth for the bar. 
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NGC 5101 [(R)SB(r)0/a] NGC5101 [(R)SB(r)0/a] 
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Fig. 1. Radial distribution of tw6 ( in c dices for each of the galaxies (a combined r ^e c fficlex and the index or the image depending 
on the emission in the bin, see text). The open circles and the filled triangles represent both sides of the bar. We have marked 
all the points with d etected emission with a star, where the detection threshold has been calculated with the prescriptions given 
in ISarzi et al.l d2006l) . When the emission completely fills all the absorption Balmer lines, the correction is more uncertain. For 
thoses cases, we have marked the points with an arrow. The vertical dashed lines represent the end-point of morphological features, 
meaning the inner bar or nuclear disk and the primary bar. Notice the general change in the indices for the two regions and the 
gradients in both. 



